Forward Neutral Pion Production in p+p and d+Au Collisions at y/s^N = 200 GeV 
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Measurements of the production of forward tt" mesons from p+p and d+Au collisions at y^sjvjv = 
200 GeV are reported. The p+p yield generally agrees with next-to-leading order perturbative QCD 
calculations. The d-|-Au yield per binary collision is suppressed as 77 increases, decreasing to ~ 30% 
of the p+p yield at {77) = 4.00, well below shadowing expectations. Exploratory measurements of 
azimuthal correlations of the forward tt" with charged hadrons at ~ show a recoil peak in p+p 
that is suppressed in d+Au at low pion energy. These observations are qualitatively consistent with 
a saturation picture of the low- a; gluon structure of heavy nuclei. 
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PACS numbers: 24.85.+p,25.75.-q,13.85.Ni,13.85.Fb 



Little is known about the gluon structure of heavy 
nuclei Q. For protons, the gluon parton distribution 
function (g-PDF) is constrained at small x (fraction of 
nucleon momentum) primarily by scaling violations ob- 
served in dccp-inclastic lepton scattering (DIS) at the 
HERA collider Q- The proton DIS data are accurately 
described by evolution equations of Quantum Chromody- 
namics (QCD) that allow the determination of the g-PDF 
[3- As cc decreases, the g-PDF is found to increase from 
gluon splitting as the partons evolve. At a sufficiently 
small value of x, yet to be determined by experiment, 
the splitting is expected to become balanced by recombi- 
nation as the gluons overlap, resulting in gluon saturation 
0. At a given x, the density of gluons per unit transverse 
area is expected to be larger in nuclei than in nucleons, 
thus, nuclei provide a natural environment in which to 
search for gluon saturation. Fixed target nuclear DIS ex- 
periments are restricted in the kinematics available; they 
have determined the nuclear g-PDF only for x 0.02 

Using factorization in a perturbative QCD (pQCD) 
framework, PDFs and fragmentation functions (FFs) 
measured in electromagnetic reactions are used to cal- 
culate hadronic processes. In p-l-p collisions, Next-to- 
Leading Order (NLO) pQCD calculations quantitatively 
describe inclusive tt*^ production over a broad range of 
pseudorapidity {i] — ~ lir[tan(6'/2)]) at ccnter-of-mass en- 
ergy = 200 GeV [10, but not at lower ^/i 0. In 
pQCD. hadroproduction at large rj from p+p collisions at 
— 200 GcV probes gluons in one proton using the va- 
lence quarks of the other, covering a broad distribution of 
gluon X peaked around 0.02 Analogously, hadropro- 
duction in the d-beam (forward) direction of d+Au col- 
lisions is sensitive to the gluon structure of the Au nu- 
cleus. Quantifying if saturation occurs at RHIC energies 
is important because the matter created in heavy-ion col- 
lisions comes predominantly from the collisions of low- 
X gluons 0. Recently, the yield of forward negatively 
charged hadrons (ft.~) in d+Au collisions was found to 
be suppressed relative to p+p 0|. The suppression is 
especially significant since isospin effects should reduce 
hr production in p-l-p collisions, but not in d+Au 1^. 

Many models try to describe forward hadroproduction 
from heavy nuclei. In the Color Glass Condensate (CGC) 
formulation, the low- a; gluon density is saturated, result- 
ing in dense color fields that scatter the partons from 
the deuteron beam 0. The average gluon- a; decreases 
rapidly with increasing 77 to « 10~^ for pions produced 
at 77 = 4 Another approach scatters quarks coher- 
ently from multiple nucleons, leading to an effective shift 
in gluon- 2; |l3l | . Shadowing models modify the nuclear g- 
PDF in a standard factorization framewo rk Other 
models include limiting fragmentation parton re- 

combination [1^ , and factorization breaking |l7j . 



Additional insight into the particle production mech- 
anism can be gained by analyzing the azimuthal corre- 
lations (A0) of the forward 7r° with coincident hadrons. 
Assuming collinear elastic parton (2 2) scattering, a 
back-to-back peak at A0 = tt is expected, with the ra- 
pidity of the recoil particle correlated with x of the struck 
gluon. In a saturation picture, the quark undergoes mul- 
tiple interactions through the dense gluon field, resul ting 
in multiple recoil partons instead of a single one [l^ Il8l |. 
thereby modifying the Ac/) distribution and possibly lead- 
ing to the appearance of monojets [l9| . 

In this Letter, we present the yields of high energy 
7r° mesons (25 < E^^ < 55 GeV) at forward rapidities 
(3.0 < 77 < 4.2) from p+p (Fig. [TJ and d+Au (Fig. ^ 
collisions at y/sJjN = 200 GeV. The data are compared 
with models and with h~ data at smaller rj. The A(j) 
distributions of the forward tt" with midrapidity are 
presented. 

Data were collected by the STAR experiment (Solenoid 
Tracker at RHIC) at the Brookhaven National Labora- 
tory Relativistic Heavy Ion Collider (RHIC). At midra- 
pidity, a time projection chamber is used to detect 
charged particles, while a forward 7r° detector (FPD) is 
used at forward rapidities. In 2002, p+P collisions were 
studied with a prototype FPD (PFPD) In 2003, p+p 
collisions were studied with the complete FPD and ex- 
ploratory measurements were made for d+Au collisions. 

The luminosity was determined using the rate of coin- 
cidences on either side of the interaction region between 
beam-beam counters (BBC) for p+p collisio ns Isll and 
zero-degree calorimeters for d+Au collisions |20j. For 
p+p, the transverse size of the colliding beams and the 
number of colliding ions were measured, giving a coinci- 
dence cross section of 26.1 ± 0.2 (stat) ± 1.8 (sys) mb |2]| . 
For d+Au, the coincidence cross section was measured to 
be (19.2 + 1.3)% of the hadronic cross section, aJ^^J^ 
The integrated luminosity for these data was w 350 nb~^ 
(200 ^b-i) for p+p (d+Au) collisions. 

Events required more energy in the calorimeter than 
from a 15-GeV electron. A BBC coincidence reduces 
non-collision background but requires an .^^-independent 
10% correction to the yields [a| to account for its effi- 
ciency. The energy is calibrated to w 1% from the cen- 
troid of the tt^ peak in the diphoton invariant mass, M^^ 
[2^ . Monte Carlo simulations with physics backgrounds 
and the full detector response describe p+p and d+Au 
data for many variables, e.g., Mjj in Fig. |21 (inset). Jet 
background is reduced in the FPD by requiring two re- 
constructed photons {Ny =^ 2), selecting 78% (53%) of 
events with > 25 GeV and iV.^ > 2 in p+p (d+Au) 
data. The 7r° detection efficiency is determined in a ma- 
trix of E'tt and rj from background-corrected simulations. 
For d+Au it is dominated by the FPD geometrical ac- 



4 




FIG. 1: Inclusive n cross section for p+p collisions versus 
the leading tt" energy (-E^) averaged over 5 GeV bins at fixed 
pseudorapidity (ry). The error bars combine statistical and 
point-to-point systematic errors. The curves are NLO pQCD 
calculations using two sets of fragmentation functions (FF). 



ceptance and is within 8-19% of the efficiency in p-l-p. 

Inclusive 7r° cross sections for p-l-p collisions at -y/s = 
200 GeV are seen in Fig.[T]at (t]) = 3.3, 3.8 0, and 4.00. 
Data are in 5 GeV bins, plotted at the average iJ^. Data 
at (77) = 3.3 and 3.8 were taken with the PFPD, where 
the systematic error increases with E^^ from 10 — 26%, 
dominated by the correction for the jet accompanying 
the 7r° 0. Data at (ry) = 4.00 were taken with the FPD, 
where the systematic error is 8 — 16%, dominated by the 
energy calibration 0. The normahzation error is 17% 
for both p+p and d-|- Au, dominated by the absolute j] un- 
certainty [221 . The curves are NLO pQCD calculations 
m using CTEQ6M PDFs and equal renormaliza- 
tion and factorization scales of pt = E^^/ cosh 77. Scale 
dependence is comparable at 77 « 4 and 77 « 0. Theoret- 
ical systematic errors, attributed to scale dependence at 

« 0, may require further study at large rj. The solid 
and dashed curves use Kniehl-Kramer-Potter (KKP) [2^ 
and Kretzer |26l | FFs, respectively, which differ primar- 
ily in the gluon-to-pion FF. Differences between FFs may 
occur at pT ^ 2 GeV/c, where the dominant contribution 
to 7r° production becomes gg scattering [23|. At (rj) = 3.3 
and 3.8, the data are consistent with KKP. At (77) — 4.00, 
the data drop below KKP and approach Kretzer as pr 
decreases, similar to the trend seen at 77 « Q. 

The study of effects from possible gluon saturation in 
a nucleus begins with the inclusive 7r° cross section for 
d-|-Au collisions (Fig.|21l. No explicit constraint is placed 
on the centrality of the collisions analyzed. The system- 
atic error is 10—22%, dominated by the background 
correction. The solid (dashed) curve is a NLO pQCD 
calculation using Au PDFs with shadowing 0| and KKP 



FIG. 2: Inclusive tt cross section per binary collision for 
d+Au collisions, as in Fig. The curves are calculations de- 
scribed in the text. (Inset) Diphoton invariant mass spectrum 
for data (stars), normalized to simulation (histogram). 

(Kretzer) FFs. The dotted curve is a LO calculation of 
multiple parton scattering [l^ , normalized to tt*' data at 
?7 « . The dot-dash curve is a LO calculation convo- 
luting CTEQ5 PDFs and KKP FFs, replacing the hard 
partonic scattering with a dipole-nucleus cross section to 
model parton scattering from a CGC in the nucleus ^3, 
normalized to d -I- Au h~ + X data at = 3.2 lOj. 
The CGC calculation overpredicts the ir'^ data here by 
a factor of 2, a factor that could approach unity with 
use of the Kretzer FF. The px dependence of the yield is 
consistent with the CGC calculation. 

The nuclear modification factor is defined as: 

oY <ei E cPo/dp^ {d + Au^Y + X) 

(^binXd" Ed^a/dp^p + p^Y + X) ■ ' 

The inelastic p-l-p cross section is af^^^ = 42 mb, while 
^hadr = (2.21 ± 0.09) b and the mean number of binary 
collisions, (i Vhin ) = 7.5 ± 0.4, are from a Glauber model 
calculation [20|. The prefactor in JS^^u equal to the 
ratio of binary collisions in p-l-p and d-|-Au, 1/(2 x 197). 
Fig. El shows i?^ A „ versus pt at (77) = 4.00 with h~ data 
at smaller 77 [lOj . Systematic errors from p-l-p and d+ Au 
are added in quadrature. The normalization error in- 
cludes the (iVbin) error but not the absolute 77 error, since 
the FPD position was the same for d+Au and p+p data. 

In the absence of nuclear effects, hard processes scale 
with the number of binary collisions and R^j^^ = 1. At 
midrapidity, R !^\.. ~ 1, with a Cronin enhancement for 
PT ~ 2 GeV / c E(l| . As 77 increases, R^j^^ becomes 
much less than 1. This decrease with 77 is qualitatively 
consistent with models that suppress the nuclear gluon 
density in III 111 Hi- Scaling i?^ by 2/3 to account 
for isospin effects on p+p — > h~+X Q , iJJ^u consistent 
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FIG. 3: Nuclear modification factor (i?dAu) for minimum- 
bias d-|-Au collisions versus transverse momentum (pt)- The 
solid circles are for n" mesons. The open circles and boxes 
are for negative hadrons fioll . The error bars are statistical, 
while the shaded boxes are point-to-point systematic errors. 
(Inset) i?dAu for tv^' mesons with the ratio of curves in Figs.|21 
andlH 

virith a linear extrapolation of the scaled i?dAu to 77 = 4. 
The curves in the inset are ratios of the calculations in 
Figs. |5] and The data lie below all the predictions. 

Exploratory measurements of the azimuthal correla- 
tions between the forward 7r° and midrapidity are 
seen in Fig. 21 for p+p and d-|-Au collisions. The lead- 
ing charged particle (LCP) analysis picks the track at 
\rifi\ < 0.75 with the highest pr > 0.5GeV/c, and com- 
putes A(j> ~ 0^0 — (pLCP for each event. The Acj) dis- 
tributions are normalized by the number of tt*^ seen at 
(77) = 4.00. Correlations near A0 = arc not expected 
due to the rj separation between the 7r° and the LCP. 
The data are fit to a constant plus a Gaussian for the 
back-to-back peak centered at A<j> = it. The fit parame- 
ters are correlated, and their errors are from the full error 
matrix. The values do not depend on . The area S un- 
der the back-to-back peak is the probability that a LCP 
is correlated with a forward tt*^. The area B under the 
constant represents the underlying event. The total coin- 
cidence probability per trigger tt^ is S + B fa 0.62 (0.90) 
for p-fp (d-|-Au), and is constant with E^^. The ratio 
S/B for p-fp docs not depend on midrapidity track mul- 
tiplicity. The peak width has contributions from trans- 
verse momentum in hadronization and from momentum 
imbalance between the scattered partons. 

A PYTHIA simulation [2^ including detector resolu- 
tion and efficiencies predicts most features of the p+p 
data PYTHIA expects S « 0.12 and B « 0.46, 

with the back-to-back peak arising from 2^2 scatter- 
ing, resulting in forward and midrapidity partons that 
fragment into the tt" and LCP, respectively The width 



FIG. 4: Coincidence probability versus azimuthal angle dif- 
ference between the forward vr'^ and a leading charged particle 
at midrapidity with pt > 0.5GeV/c. The left (right) column 
is p-fp (d-|-Au) data. The curves are fits described in the text, 
including the area of the back-to-back peak (S). 



of the peak is smaller in PYTHIA than in the data, which 
may be in part because the predicted momentum imbal- 
ance between the partons is too small, as was seen for 
back-to-back jets at the Tevatron [soj . 

The back-to-back peak is significantly smaller in d-f Au 
collisions than in p+p, qualitatively consistent with the 
monojet picture arising in the coherent scattering [l^ 
and CGC [H models. HIJING [HI includes a model of 
shadowing for nuclear PDFs. It predicts that the back-to- 
back peak in d+Au collisions should be similar to p+p, 
with S w 0.08. The data are not consistent with the 
HIJING expectation at low E^,-. 

In conclusion, the inclusive yields of forward 7r° mesons 
from p+p collisions at -^/s = 200 GeV generally agree 
with NLO pQCD calculations. However, by (77) = 4.00, 
the spectrum is found to be harder than NLO pQCD, 
becoming suppressed with decreasing px- In d+Au col- 
lisions, the yield per binary collision is suppressed with 
increasing 77, decreasing to ~ 30% of the p+p yield at 
(77) — 4.00, well below shadowing and multiple scatter- 
ing expectations, as well as exhibiting isospin effects at 
these kinematics. The pT dependence of the d+Au yield 
is consistent with a model which treats the Au nucleus as 
a CGC. Exploratory measurements of azimuthal correla- 
tions of the forward tt'^ with charged hadrons at midra- 
pidity show a recoil peak in p+p collisions that is sup- 
pressed in d+Au at low Ej^, as would be expected for 
monojet production. These effects are qualitatively con- 
sistent with a gluon saturation picture of the Au nucleus, 
but cannot definitively rule out other interpretations. A 
systematic program of measurements, including direct 
photons and di-hadron correlations over a broad range of 
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A?7, pt, and ^/s, is needed to explore the nuclear modifi- 
cations to particle production. A quantitative theoretical 
understanding of the observables is needed to facilitate 
experimental tests of a possible color glass condensate. 
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